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ABSTRACT 

In this paper, we analyze in detail with numerical simulations how the bad data removal, i.e., the mask 
effect, can influence the peak statistics of the weak lensing convergence field reconstructed from the shear 
measurement of background galaxies. It is found that high peak fractions are systematically enhanced due to 
the presence of masks, the larger the masked area, the higher the enhancement. In the case that the total masked 
area is about 13% of the survey area, the fraction of peaks with signal-to-noise ratio v>3is~ 13% of the total 
number of peaks, in comparison with ~ 9% of the mask-free case in our considered cosmological model. This 
can have significant effects on cosmological studies with weak lensing convergence peak statistics, inducing 
a large bias in the parameter constraints if the effects are not taken into account properly. Even for a survey 
area of 9 deg , the bias in (D.m, os) is already intolerably large and close to 3cr. It is noted that most of the 
affected peaks are close to the masked regions. Therefore excluding peaks in those regions in the peak statistics 
can reduce the bias effect but at the expense of loosing usable survey areas. Further investigations find that 
the enhancement of the number of high peaks around the masked regions can be largely attributed to the fewer 
number of galaxies usable in the weak lensing convergence reconstruction leading to higher noise than that of 
the areas away from the masks. Based on Fan et al. (2010), we thus develop a two-noise-level model that treats 
the areas close to and away from the masked regions separately. It is shown that the model can account for 
the mask effect on peak statistics very well, and the bias in cosmological parameters is significantly reduced if 
this model is applied in the parameter fitting. We also analyze the systematic effects on peak statistics resulting 
from the nonlinear convergence reconstruction, including the smoothing order problem and the mass-sheet 
degeneracy. 

Subject headings: cosmology - dark matter - clusters: general - gravitational lensing: weak - large-scale struc- 
ture of universe 



1. INTRODUCTION 

Gravitationally induced weak lensing effects have emerged 
as one of the most important probes in cosmological studies 
(e.g., Bartelmann & Schneider 2001; Albrecht et al. 2006; 
Amendola et al. 2012; Abate et al. 2012; Heymans et al. 
2012; Erben et al. 2012; Simpson et al. 2012; Kilbinger et 
al. 2012). Besides the shear two-point correlation analyses, 
weak lensing peak statistics can provide important and com- 
plementary information especially considering that the struc- 
ture formation is a nonlinear process (e.g.. White et al. 2002; 
Hamana et al. 2004; Tang & Fan 2005; Hennawi & Spergel 
2005; Marian et al. 2009; Dietrich & Hartlap 2010; Kra- 
tochvil et al. 2010; Marian et al. 2012; Hilbert et al. 2012). 
Current generation of observations have proved the feasibility 
of detecting massive clusters from weak lensing peak identi- 
fications (e.g., Wittman et al. 2006; Gavazzi & Soucail 2007; 
Shan et al. 2012). Future weak lensing surveys will be able 
to provide a large number of peaks with high signal-to-noise 
ratio and therefore their statistics should expectedly be able to 
contribute significantly to precision cosmological studies. On 
the other hand, it is known that many effects can affect pro- 
foundly the weak lensing peak statistics. The complex mass 
distribution of clusters of galaxies and the projection effects of 
large-scale structures along line of sights prevent us from link- 
ing weak lensing peaks to single clusters in a simple way (e.g.. 
Tang & Fan 2005; Yang et al. 201 1; Yang et al. 2012; Hamana 
et al. 2012). The intrinsic ellipticities of source galaxies gen- 
erate large noise that not only produces false peaks through 
their chance alignments (e.g.. Van Waerbeke 2000; Fan 2007) 
but also affects the true peak signals from massive clusters 



significantly (Fan et al. 2010). Furthermore, various obser- 
vational effects can also have large impacts on weak lensing 
peak statistics if they are not taken into account properly. The 
full realization of the power of weak lensing analyses in future 
cosmological studies relies on our thorough understandings 
about different systematics. 

Weak lensing observations target at far way background 
galaxies, and bad data occurrences are unavoidable and they 
should be masked out carefully (e.g., Heymans et al. 2012; 
Erben et al. 2012). These masks can occupy ~ 10% to 
~ 20% of the total survey area and result irregular survey 
boundaries and artificial voids in the background galaxy dis- 
tribution, which in turn can affect the weak lensing analy- 
ses considerably. In Hikage et al. (2011), they investigate 
the mask effect on the shear power spectrum estimation and 
propose the pseudo-spectrum method to account for the ef- 
fect. Considering convergence peaks, because the conver- 
gence field is reconstructed from the shape measurements of 
background galaxies, the devoidness of galaxies in masked 
regions affects the convergence calculation and consequently 
the peak statistics inevitably. In this paper, we perform de- 
tailed studies about the mask effect on weak lensing conver- 
gence peak statistics and further the derived cosmological pa- 
rameter constraints. We adopt the simulations from White & 
Vale (2004) as the fiducial inputs for weak lensing signal cal- 
culations. Background galaxies with intrinsic ellipticities are 
randomly populated and 'observed' ellipticities including the 
shear signals from the fiducial inputs are then constructed for 
each galaxies. The masks are generated according to the mask 
size distribution from Shan et al. (2012), and are given spa- 
tial positions randomly in our statistical analyses. We then 



remove galaxies inside the masks and reconstruct the conver- 
gence field from remaining galaxies. The peak statistics is an- 
alyzed and compared with the case without masks. We further 
explore different methods to correct the mask effects. Along 
the way, we also analyze the systematic effects arising from 
the nonlinear convergence reconstruction on peak statistics. 

The rest of the paper is organized as follows. In §2, we in- 
troduce the basics related to our studies, including the conver- 
gence reconstruction method and the theoretical peak model 
of Fan et al. (2010) used in our analyses. In §3, we describe 
the different simulations constructed for our analyses and the 
procedures for cosmological parameter fitting assuming the 
peak counts from simulations as the 'observational data' . The 
main results are presented in §4. Discussions are given in §5. 

2. THEORETICAL BASICS 
2.1. Convergence reconstruction 

Observationally, the weak lensing effect is mostly extracted 
from the shape distortion measurements of background galax- 
ies, which is directly related to the weak lensing shear com- 
ponents. On the other hand, the lensing convergence is the 
weighted projection of the density distribution along the line 
of sight, and thus more visually linked to the structures we are 
interested in. In this paper, we concentrate on the convergence 
peak analyses. 

Considering small source galaxies, their linear-order im- 
age distortion from the gravitational lensing effect of a single 
lens can be described by the Jacobian matrix given by (e.g., 
Schneider et al. 1992) 
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where k is the lensing convergence and yi and 72 are the two 
shear components with 

The lensing potential ifr is determined by the surface mass den- 
sity of the lens through 
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where S^r is the critical surface mass density given by 
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with Di, Ds and D;, being the angular diameter distances from 
the observer to the lens, to the source, and between the lens 
and the source. It can be seen that k - I./'Zcr- For weak lens- 
ing effects from large-scale structures beyond a single lens, 
under the Born approximation the above formulations still 
hold except the lensing convergence is given by, in the case of 
a fixed source position (e.g., Bartelmann & Schneider 2001), 
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where w is the comoving radial distance, fj^ is the angular 
diameter distance, a is the scale factor of the universe, and 6 
is the density perturbation along the line of sight. 

The image distortion is then described by k and 7, with the 
quantity (detA)"' - [(1 - /c)^ - I7P]"' giving rise to the flux 



magnification (I7I = (y^ + 72)'^^), and the eigen values of 
A related to the axial length. Specifically, the lensing effect 
makes a circular source appear as an ellipse with the axial 
ratio of 

a^ ^ l-K-\r\ ^ 1-lgl 

b^ l-K+M l+\g\' 

where gi = 7,7(1 - k) is named as the reduced shear compo- 
nent. Thus for ideally circular sources, we can obtain g by 
accurately measuring the shape of the sources, and further re- 
construct the convergence k from the relation between k and 
7, which in the Fourier space can be written as (e.g.. Kaiser & 
Squires 1993) 
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where D is given by 
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However, galaxies have intrinsic ellipticities. The complex 
ellipticity of the lensing distorted image e with |e| = (1 - 
b/a)/(l + bid) is then related to the intrinsic ones e, by the 
following relation (e.g., Seitz & Schneider 1997) 
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for \g\ > 1 



It has been shown that the average of c gives rise to the unbi- 
ased estimate of g and l/gfoi\g\ < 1 and|g| > 1 , respectively 
(Seitz & Schnieder 1997). We further construct the distortion 
6hy (e.g., Schneider & Seitz 1995) 
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as the observed quantity, which is independent of |g| < 1 or 
> 1. One can then solve for 7 by 
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where the sign is determined by -sign[det(A)]. We proceed 
the reconstruction of the lensing convergence iteratively from 
the following relation 
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where 0(0) = (O] - Gj 



210^62) /\0\'^. Specifically, we start 
by assuming k^*^^ = and |g-| < 1 everywhere, and thus (e.g., 
Bartelmann 1995) 
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At n-th step, we obtain /c*^"* from y^" '' via Eq.(12) and further 
calculate 7,"', from a-*"' to determine the sign of det(A*"') ev- 
erywhere. At n + 1 step, we insert k^"^ into Eq. (1 1) to estimate 
-y(n) by considering the signs of det(A*"') calculated in step n. 
In the case of a- << 1 and \y\ « 1, we have (e) = -g ~ -7, 
and thus the convergence reconstruction is a single-step linear 
process. 



2.2. Peak model 

To constrain cosmological parameters from weak lensing 
peak abundances, we need to calculate the expected peak 
numbers for different cosmological models. From previous 
section, we see that the intrinsic ellipticities of lensed galax- 
ies are the dominant source of errors in weak lensing analyses. 
With the process of averaging over a number of galaxies to ob- 
tain the estimate of the lensing signal, the residual noise is on 

the order of cr^l JngO^, where cr^ is the rms of the intrinsic 

ellipticities, Ug is the surface number of lensed galaxies in the 
weak lensing analyses, and Oq is the smoothing scale used in 
the averaging. Such a noise can affect significantly the weak 
lensing peak statistics. In our studies here, we employ the the- 
oretical model of Fan et al. (2010) to calculate peak numbers 
for relatively high peaks in cosmological parameter fittings. 
The model takes into account the effects of noise by consid- 
ering the noise-induced bias and the dispersion on the height 
of the true convergence peaks from massive halos and the en- 
hancement of the noise peak abundances due to the existence 
of the true mass distribution properly. 

Specifically, the model assumes that the derived conver- 
gence field after a smoothing can be written as K^ = K + N, 
where K represents the true lensing convergence, and A^ is for 
the residual noise from intrinsic ellipticities. The noise field A^ 
is modeled as a Gaussian random field from the central limit 
theorem (e.g.. Van Waerbeke 2000). The model concentrates 
on high peaks and assumes that true peaks come from indi- 
vidual massive halos. Thus a considered area is split into halo 
regions and field regions. In an individual halo region, the 
peak number distribution can be calculated from the Gaussian 
statistics of Kfj with known K from the halo. Then the to- 
tal number of peaks in halo regions can be obtained from the 
summation of the peaks in halo regions weighted by the halo 
mass function. In the field region, the numbers of peaks are 
directly computed from the noise field A^. The total surface 
number density of peaks is then written as 



nj,eakiv)dv = rfp,„u(v)dv + nl^^^{v)dv. 
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where v = Kn/ctq is the signal-to-noise ratio of a peak and 
the term «'^g„j.(v) is for halo regions including both true peaks 
corresponding to real halos and the noise peaks within the 
halo regions, and n" ^^^iy) is for field regions with only noise 
peaks. For «'^g^^(v), We have 
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where dV(z) is the cosmological volume element at redshift z, 
dO. is the solid angle element, «(M, z) is the mass function of 
dark matter halos and. 
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gives rise to the number of peaks in the area within the virial 
radius of a halo of mass M at redshift z- Here Upeakiv, R, M, z) 
describes the surface number density of peaks at the location 
of R from the center of the halo, which depends on the pro- 
jected density profile of dark matter halos and is given by 
Eq.(15) of Fan et al. (2010). We adopt the Sheth-Tormen 
mass function in our calculations ( |Sheth & Tormen"||1999| l. 
The linear power spectrum used in computing the rms den- 
sity perturbations is the same as that used in the simulations 



of [White & Wel ( p004l l from jEisenstein & Hu | ( fT999| . The 
NFW profile (Navarro et al. 1996, 1997) is appUed to de- 
scribe the density distribution of dark matter halos with the 
concentration-mass relation given by Bullock et al. | ( [200T] ) 

c„>(M,z)=-^ — — — - , (17) 
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where c» is taken to be c» = 8. The field term ^"^.^^(v) is given 
by 
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where n,an{y) is the surface number density of pure noise 
peaks (Fan et al. 2010). 

3. SIMULATIONS AND COSMOLOGICAL 
PARAMETER FITTING FROM PEAK WEAK 
LENSING PEAK COUNTS 

3.1. Base simulations 

In our analyses, we use a set of lensing maps constructed 
from a flat ACDM N-body simulation by [White & Vale 
(2004 ) as our base inputs. The cosmological parameters of the 
considered simulation are Q.„ = 0.296, h = 0.7, cr^ = 0.93, 
and n - \, where Q,„ and h are the present dimensionless mat- 
ter density of the universe and the Hubble constant in units of 
lOOkms^^Mpc"', respectively, erg is the rms of the extrap- 
olated linear density perturbations smoothed over 8/i"'Mpc, 
and n is the power index of the initial power spectrum of den- 
sity fluctuations. The N-body simulation box is 200/i"'Mpc, 
and the force softening length is ~ 30/; 'kpc. The simula- 
tion starts from z - 50, and the initial conditions are gener- 
ated from the transfer function for the power spectrum given 
by Eisenstein & Hu (1999) with the baryonic oscillations 
smoothed out (White & Vale ,200 4). Two sets of simula- 
tions are done with different reaUzations of the initial con- 
ditions. The lensing convergence and shear are calculated 
through the ray-tracing technique. The source redshift is set 
to be Zs = 1. With 16 from each simulation, we have total 
32 sets of 3 X 3 deg" converge nce and she ar maps constructed 
along different line of sights ( [White &"V ale 2()04), and each 
has 1024 X 1024 pixels. They form the base data sets for our 
weak lensing peak studies. It is noted that th e cosmological 
param eters are in accord with WMAPl (e.g., Spergel et al. 
2003 | l and different somewhat from the most rece nt fitting re- 
sults, es pecially for (Til, (e.g, Komatsu et al. |20TT|[Hinshaw et 



[al. |2012) . However, because we mainly focus on analyzing 
mask effects by comparing weak lensing peak statistics with 
and without masks, the choice of cosmological model should 
not affect our results significantly. 

3.2. Fiducial reconstructed convergence maps 

To analyze weak lensing convergence peaks, previous stud- 
ies often add a Gaussian noise arising from intrinsic elliptici- 
ties directly to convergence map s from N-body sim ulations to 
construct 'observed' maps (e.g.,[Hamana et al. |20 04a; Tang] 
& Fan |[2005| [Yang et al. l2012[ ). ObservationaUy however, 
a convergence map has to be reconstructed from the mea- 
sured galaxy ellipticities. To mimic observations, we thus also 
perform the convergence reconstruction. We generate source 
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Figure 1. Upper-Left: One example of the smoothed base convergence map without noise; Upper-Right: The corresponding shear reconstructed convergence 
map with intrinsic elHpticities of source galaxies included; Lower-Left: The corresponding g reconstructed convergence map; Lower-Right: The corresponding 
masked g reconstructed convergence map, here the yellows are the masks in this case. 



galaxy catalogs by randomly populating them in angular po- 
sitions at Zs = 1 and assigning them intrinsic elHpticities ac- 
cording to the following probability distribution 



Ps{^s\,e,2) 
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TTo-Jl -exp(- 1/0-2)]^ 



\es\ e [0, 1] (19) 



where e^i and esi are the two components of the intrinsic elHp- 
ticities, le,! - Je^^ + e^2 ^^^ ^h^ ^m^ dispersion of le^l is taken 
to be cTf: - 0.4. We assume the number density of source 
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galaxies to be «g = 30 arcmin" , and the spatial clustering 
and the intrinsic alignment of source galaxies are not consid- 
ered here. We calculate the reduced shear signal g for each 
source galaxy from the base pure lensing maps by interpola- 
tion from regular grids to galaxy positions. The 'observed' 
galaxy elHpticities e are then constructed according to Eq.(|9]l. 
We calculate the average (e) on positions of regular grids as 



follows 
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where the summation is over galaxy positions represented by 
the subscript /. The window function W is taken to be Gaus- 
sian given by 
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The smoothing scale 0a is set to be 6*0 = 1 arcmin. With 
the obtained (e), the convergence reconstruction is done itera- 
tively as described in §2.1. The results converge quickly with 
about eight iterations for the converging accuracy of 10"''. For 
each base map, we generate 5 sets of galaxy catalogs with 
different realizations of spatial positions and intrinsic elHp- 
ticities. We thus obtain total of 32 x 5 = 160 reconstructed 
3x3 deg^ convergence maps. We refer such maps as 'g re- 
construction' . For comparison, we also generate sets of 'ob- 



served' ellipticities by e = ^s - 7 and perform linear recon- 
struction to get convergence maps, which are referred to as 
'shear reconstruction' maps. 

Figure 1 shows a set of convergence maps. The upper left 
is the base pure convergence map smoothed with a Gaussian 
window function with Qq - \ arcmin. The upper right and 
lower left are the reconstructed maps from shear and from the 
reduce shear g, respectively. We can see that the reconstructed 
maps are noisy comparing to the base one, but most of the 
high peaks are still apparent. On the other hand, the noise 
affects the peak height and the pure noise peak distribution 
is also biased by the true mass distribution, which have to 
be taken into a ccount properly in modeling the weak lensing 
peak statistics ( [Fan et al. [2010 1. The two reconstructed maps 
from shear and from g are very similar, but there are subtle 
systematic differences, for which we will discuss later in §4. 
The lower right map is the reconstruction with masks to be 
described in the next subsection. 

3.3. Mask model and convergence reconstruction with 
masks 

Removing bad/low quality imaging data is essential in weak 
lensing observational analyses. This leaves holes in the source 
galaxy distribution, which in turn affects the convergence re- 
construction and subsequent cosmological studies. To inves- 
tigate the mask effects on weak lensing peak counts statisti- 
cally, we generate mock masks by modeling the basic masks 
for point sources, bright saturated stars and bad pixels being 
circular in shape. The mask si ze distribution is in accord with 
that of CFHTLS used in Sha n et al. ] p012] l. We also add 
rectangle-shaped masks in both x and y directions to those 
circular ones with radius larger than 1 arcmin to mask out sat- 
uration spikes. These extra masks have a size of 0.2rx5r with 
r being the radius of the circular mask to be added on. We 
populate masks randomly in each of the considered 3x3 deg 
fields. With the size distribution of |Shan et al. | ( |2012| ), the 
total number of masks in each field is set to be Nmask- We con- 
sider three cases with N,„ask - 140, 280, and 420, correspond- 
ing to the total masked area fraction of ~ 6%, ~ 13% and 
~ 19%, respectively. We then remove galaxies within masks 
from the source galaxy catalogs generated in §3.2 and recon- 
struct convergence maps from the remaining source galaxies. 
An example is shown in the lower right panel of Figure 1 . The 
masks are shown in yellow. We then have two separate sets of 
convergence maps reconstructed from 'observed' ellipticities 
without and with masks, respectively, and each set contains 
totally 160 of 3 x 3 deg^ convergence maps for peak analyses. 

3.4. Cosmological parameter constraints 

The mask effects on cosmological studies with weak lens- 
ing peak statistics are investigated by analyzing their impacts 
on cosmological parameter constraints from weak lensing 
peak counts. As described in §2.2, |Fan et al7~| (2010) pro- 
pose a model to calculate the peak counts including both the 
true peaks from massive halos and noise peaks from chance 
alignment of background galaxies by taking into account the 
influence of noise on true peaks and the noise peak enhance- 
ment due to the mass distribution around halos properly. The 
model allows us to use directly the number of peaks detected 
from 'observed' convergence maps for cosmological studies 
without the need to identify true peaks associated with mas- 
sive halos. It should be pointed out that the model of Fan et 
al. (2010) is built on the assumption that true peaks are asso- 
ciated with individual dark matter halos. It has been shown 



that t his is largely valid for high peaks (e.g., Yang et al. 
2011 1. For low peaks, on the other hand, the projection ef 
feet from large scale structures along line of sights is sig- 
nificant. We are in the process developing a more complete 
model including the projection effects. With that, the cos- 
mological constraints from peak counts can be extended to 
include low peaks to further enhance the cosmological gains. 
At the moment, however, we only concentrate on high peaks. 
For the noise level considered here, we include peaks with 
3.75 < V < 6.75 in the analyses, where the signal-to-noise 
ratio of peaks is defined as v = Kf^/cro with ctq the rms of 
noise from intrinsic ellipticities. In t he Gaussian smoothing 
used here, crii - crJ(47Tng6^y'^ (e.g. van Waerbeke [20001. 



For o-f: - 0.4, ng = 30 arcmin" and 6c - i arcmin, we have 
cr^) = 0.02. The high-end cut is introduced to avoid the possi- 
ble bias from the very limited number of highly massive halos 
existed in the relatively small-sized simulations we used. 

We construct the 'observed data' specifically as follows. 
For each of the considered case, we have 160 of 3 x 3 deg^ 
convergence maps. To reduce the boundary effects on peak 
counts, we exclude the outermost 10 pixels in each direc- 
tion in peak counting. Thus the effective area of each map 
is [3(1 - 20/1023)]^ ~ 8.65 deg^. For each map, we count 
peaks in 6 signal-to-noise ratio bins of width 0.5 in the range 
of 3.75 < V < 6.75. We further obtain the mean number of 
peaks in each bin by averaging over the 160 maps. We then 
scale the peak counts back to 9 deg" by multiplying a factor 
of 9/8.65 on the peak counts in each bins. These peak counts 
form the 'observed' data denoted by A^,, (i - 1, ..., 6). For cos- 
mological parameter fitting, we also calculate the covariance 
matrix by 

^'■^- = ^TT TS^i" - ^^^^T - ^>^' ^22) 

where r denotes different maps and R - 160. The parameter 
constraints are then obtained by minimizing the x^ defined by 

XJ, = dN^^'\cr^dN^"'> = Yj dNf\C;/)dNf\ (23) 

0=1,...,6 

where dN- - Nf -Ni with N^ being the model prediction 
for the cosmological model /' . 

4. RESULTS 

As described previously, we use the model of Fan et al. 
(2010) as the basic theoretical model. Thus before presenting 
the mask effects, we first need to show the cosmological ap- 
plicability of the model. Although it has been de monstrated 
that th e model works well for high peak counts ( Fan et aiT~| 
2010| l, here we test the model applicability quantitatively by 
performing cosmological parameter constraints from the 'ob- 
served' peak counts without masks. Along the way, we also 
discuss the effects of reconstruction on peak counts. With 
these analyses results, we then proceed to discuss the mask 
effects. 

4.1. The model applicability in cosmological parameter 
constraints 

To exam the model applicability and possible effects aris- 
ing from convergence reconstruction, we construct separate 
sets of 'observed data' from convergence maps obtained by 
different ways. In the very weak lensing regime, the observed 
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Figure 2. Upper-Left: The average peak counts in 1 deg". The blue, red and black histograms are for the results from direct noise-added convergence maps, 
from the model prediction, and from the noise-free base convergence maps. Upper-Right: Cosmological constraints on (fi,,,, irg) from^- fitting for a survey area 
of 9 deg for the direct noise-added case. The solid lines are lo", Icr and 3cr confidence contours. The blue symbol represents the fiducial input parameters of the 
simulations and the red symbol is the best-fit result. The dashed line shows the best fit erg given n„, and represents more or less the direction of the degeneracy 
of the two parameters. Lower-Left: Same as the upper-right panel but for the survey area S = 100 deg". Lower-Right: Same as the upper-right one, but for the 
survey area S = 1 8, 000 deg~ . 



ellipticity for a background galaxy can be written as e ^ Cj-7, 
and consequently, the reconstructed convergence after a suit- 
able smoothing is K^ = K + N where K is the smoothed lens- 
ing convergence and A^ is the residual noise due to the finite- 
number averaging of the intrinsic ellipticities of source galax- 
ies, which is approximately Ga ussian due to the central limit 
theorem (e.g., [van Waerbeke |2000j . Therefore a straightfor- 
ward way to include the noise from intrinsic ellipticities into 
consideration is to linearly add the noise field A^ directly to the 
base convergence maps from simulations. This approach has 
been adopted by many previous numerical studies on weak 
lensing peak statistics (e.g., Hamana et al. 2004a Tang & 



Fan |2005[|Yang et al. 12012^ . As discussed in §2.2,' the form 
oi Kn = K + N is also the starting point of the model of Fan et 
al. (2010). We thus first construct a set of noisy convergence 
maps by adding a random Gaussian noise to the 1024 x 1024 



grids of the base convergence maps with the variance given 
by (e.g., Hamana et al. |2004a 1 



o-J 



«*e2,; 



(24) 



where o-^ 



0.4, 



30 arcmin 



and 6„ 



180/1024 arcmin is the angular size of a pixel. We then 
smooth the noise-add convergence maps with a Gaussian 
smoothing of flg = 1 arcmin. We generate 5 noisy con- 
vergence maps with different realizations of the noise field, 
and perform peak counting from the total of 160 convergence 
maps to obtain the first set of 'observed data' and the corre- 
sponding covariance matrix. Note that the model of Fan et al. 
(2010) is a spherical halo-based model with only noise effects 
included. Thus by carrying out cosmological parameter fit- 
ting to this set of 'observed data', we can test if other effects 




Figure 3. Cosmological parameter constraints from x^ fitting for 9 deg' 
shear reconstruction case. The lines and symbols are the same as those in 
the upper right panel of Figure 2. 



beyond the noise from intrinsic ellipticities, e.g., the projec- 
tion effects from large-scale structures and the non-sphericity 
of dark matter halos, can affect weak lensing peak counts and 
therefore bias the fitting result significantly. 

The fitting results are presented in Figure 2. The upper left 
panel shows the peak counts per deg^ from the noisy conver- 
gence maps (blue), the Fan et al.'s model prediction (red) and 
the peak counts from base convergence maps without noise. 
It is seen that the existence of noise systematically increases 
the number of peaks due to its influence on true peak heights 
and the number enha ncement o f pure noise peaks around dark 
matter halo regions (Fan et al. ,2010) . The model prediction 
agrees with the results from noisy convergence maps quite 
well. It is noted that the base convergence maps are gener- 
ated from two sets of N-body simulations with the box size 
of 200/i"'Mpc. Because of the relatively small simulation 
box, the number of massive dark matter halos is limited and 
the sample variance can be significant. Although we have 
excluded the very high peaks in the consideration, the peak 
counts from the two sets of 16 base convergence maps from 
two different realizations of the simulations are apparent. This 
sample variance contributes dominantly to the differences be- 
tween the blue and the red histograms at high ends. The good 
agreement between the model prediction and the 'observed 
data' shows that the noise effects from intrinsic ellipticities 
of background galaxies are the leading systematics for weak 
lensing peak counts in the considered range and the projection 
effects and the non-sphericity of dark matter halos play minor 
roles statistically. The upper right panel shows the fitting re- 
sult with the observed area of 9 deg^. The covariance matrix 
used in the fitting is calculated directly from the 160 noisy 
convergence maps. The blue symbol is the fiducial parame- 
ter input of the simulation with Q„, = 0.293 and cr-i, - 0.93, 
and the red symbol is the best fit based on the model of Fan 
et al. (2010). The solid lines from inside out are the Icr, 
2cr and 3cr contours, respectively. The dotted line shows the 
best fit cr% given Q„, and represents more or less the degener- 
acy between the two cosmological parameters with respect to 
the peak counts considered here. The location of the best fit 





Figure 4. The normalized covariance matrix in grey scale for the direct 
noise-added case (left) and for the shear reconstruction case (right). 



(fim, (Tg) along the dotted line changes somewhat from one set 
of runs to another, again owing to the degeneracy of the two 
parameters and the y^ differs only slightly along the dotted 
line. The fact that the blue symbol is very close to the dotted 
line shows that the model used in the parameter fitting works 
well without introducing significant bias. We also consider 
the cases with the survey area of 100 deg^ and 18000 deg^ 
shown in the lower left and lower right panels, respectively. 
We take a simple approach to rescale the covariance matrix 
calculated from 160 noisy convergence maps to the two larger 
surveys as suming a Poisson scaling relation to the survey area 
S as 1/5 ( Kratochvil et al. |2010| l. This may underestimate 
the covariance matrix by a factor of ~ 1.5 give n the exis- 
tence of long -range correlations of the true peaks fKratochvil 
|et al. pOlOp . Furthermore, for our analyses here, the base 
convergence maps are derived from only two sets of simula- 
tions with a relatively small simulation box. Thus the sample 
variance cannot be averaged out, which affects the 'observed' 
peak counts as shown in the upper left panel of Figure 2. This 
also leads to an underestimate of the covariance matrix. Keep- 
ing these in mind, the fitting result for S - 100 deg^ is still 
good without significant bias comparing to the statistical er- 
rors. For the case with S - 18000 deg^, the best fit values 
are away from the fiducial inputs significantly in a statisti- 
cal sense. Thus for future precision cosmological studies, we 
need to further develop our model to take into account more 
effects in addition to the noise effects. On the other hand, as 
discussed above, to fully quantify the model performance, we 
also need larger simulations so that the sample variance bias 
can be well controlled and the covariance matrix can be suit- 
ably estimated. 

The main focus of the paper is to investigate the mask ef- 
fects on weak lensing peak counts. For that purpose, we con- 
clude that the model of Fan et al. (2010) is good enough to 
reveal the bias on cosmological parameter fitting induced by 
the mask effects because of their significance (see §4.2). 

To be more close to observational analyses, we also con- 
struct two sets of 'observed data' and the corresponding co- 
variance matrix from convergence maps obtained from shear 
reconstruction and reduced-shear g reconstruction, respec- 
tively. For the shear reconstruction described in §3.2, we ar- 
tificially generate 'observed ellipticities' of source galaxies 
by e = fj - 7 and then reconstruct the convergence maps 
linearly. Mathematically, we expect that the obtained con- 
vergence maps should be exactly the same as Kj^ - K + N. 
The fitting result is shown in Figure 3. Noting the approxi- 
mate degeneracy of the two parameters along the dotted line, 
the best fit result is indeed very consistent with that shown in 
the upper right panel of Figure 2. However, the confidence 
contours here are considerably tighter. Shown in Figure 4, 
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Figure 5. The left panel shows the average peak counts before (blue) and after (red) the mass-sheet coiTection for the shear reconstniction case. The right panel 
shows the corresponding cosmological parameter constraints after the correction. The blue, red and green symbols represent the fiducial input of the simulations, 
the best fit after and before the mass-sheet connection, respectively. 
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Figure 6. The results for the g reconstruction case. The left panel shows the average peak counts with the blue one from g reconstructed maps and the red one 
for the theoretical predictions. The right panel shows the corresponding cosmological constraints. The red, green and blue symbols are the best fit from the g 
reconstruction, the shear reconstruction and the fiducial input, respectively. 



we exam the covariance matrices from the direct noise-added 
maps (left) and the shear reconstructed maps (right), and find 
that the cross correlations are much smaller in the latter case. 
Careful analyses show that this difference is due to the mass- 
sheet degeneracy inherited from the reconstruction. In the re- 
construction case, the average {K^) - over the whole map. 
For the direct noise-added case, however, {K^} x (K). For a 
region of 3 x 3 deg^, the true average K is on the order of 10"'' 
and different maps have different {K}. The fluctuations in {K} 
lead to larger cross correlations in the covariance matrix cal- 
culated from the direct noise-added maps than that in the re- 
construction case with {K^} - for all the maps. On the other 
hand, because the fluctuations in {K} can be either positive or 
negative, there is no net effect on the average peak counts. 
To demonstrate this, we make corrections to each shear re- 
constructed map by adding the value (K) of the correspond- 



ing base convergence map to it and then perform the peak 
counting and calculate the covariance matrix. The results are 
shown in Figure 5, where the left panel shows the average 
peak counts from shear reconstructed maps before (blue) and 
after correction (red), and the right panel shows the parameter 
constraints from the corrected shear reconstructed maps with 
the blue symbol for the fiducial input, the red symbol for the 
best fit from the corrected maps and the green one for the best 
bit from the uncorrected shear reconstructed maps, it is seen 
that the confidence contours are indeed enlarged and are very 
much the same as that of the upper right panel of Figure 2. It 
should be pointed out that this mass-sheet degeneracy effect 
depends on the considered survey areas, the larger the area, 
the weaker the effect. 

We now turn to the results from g reconstruction, which 
mimics the observational analyses the most. The results are 
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Figure 7. Test results for a known dark matter halo. The upper left panel shows the smoothed convergence profiles from theoretical calculations (solid blue), the 
shear reconstruction (dashed red) and the g reconstruction (dash-dotted green), respectively. The lower left panel shows the corresponding relative differences. 
The right panels are the same as the left one except the g reconstruction results are corrected. 
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Figure 8. Similar to Figure 6, but for the results of corrected g reconstruction. 

shown in Figure 6. The left panel presents the average peak 
counts from the g-reconstructed maps (blue) and the model 
prediction of Fan et al. (2010). It is seen that the blue his- 
togram is systematically lower than the red one. This is also 
reflected in the cosmological parameter fitting shown in the 
right panel. The fiducial input (blue symbol) is clear above 
the degeneracy dotted line where the best fit (red) is on al- 
though the difference is still well within the Icr range here. 
For comparison, we also show the best fit from the shear re- 
construction case (green symbol), which is closer to the fidu- 
cial input considering the degeneracy of the two parameters 
(see Figure 3). Similar to the shear reconstruction case, here 
the mass sheet degeneracy also causes apparently tighter con- 
fidence levels than those of the direct noise-added case. 

For the systematically lower peak counts from g reconstruc- 
tion, we find that it results largely from the smoothing in g. As 
we described in §3.2, to carry out the reconstruction, we first 
need to obtain a local estimate of g by averaging (smoothing) 
over surrounding source galaxies. Such a procedure gives rise 



to {g} = (7/(1 - k)} (in the sub-critical case with |g-| < 1), 
but not {y)l{\ - (k)). Consequently, the reconstructed con- 
vergence is not the same as the real smoothed convergence. 
Specifically, the real smoothed convergence should be mathe- 
matically obtained by 
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(25) 
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The difference between {g} and {y)l{\ - {k)) (denoted by a) 
leads to the difference between the g reconstructed conver- 
gence (denoted as {Kg) for the moment) and the real {k). Let 
Ax - {Kg) - {k), we then have, to the Hnear order in a and Ajc, 
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(26) 



We demonstrate this smoothing order problem by consider- 
ing the weak lensing signal from a single known NFW dark 
matter halo with mass of M - 10'^ M©. The left panels 
of Figure 7 show the results by setting the intrinsic elliptic- 
ities of source galaxies to be zero to isolate the smoothing 
order problem. The upper panel shows the smoothed conver- 
gence profiles from the theoretical calculations (solid blue), 
the shear reconstruction (dashed red) and the g reconstruction 
(dash-dotted green), respectively. The lower panel gives the 
relative differences between the theoretical result and the re- 
sults from shear reconstruction (blue) and the g reconstruction 
(dashed red), respectively. It is seen that while there is almost 
no difference for the shear reconstruction case, the relative 
difference is on the order of 1% in the central region of the 
halo in the g reconstruction case. This results the peak height 
from g reconstruction to be lower than the theoretical one by 
Q.Q\{kq). Because the dark matter halo is known, we can cal- 
culate a, and further to get Ax iteratively. The right panels 
of Figure 7 show the corrected results for g reconstruction. 
We can see that after the correction, the difference between 
the theoretical calculation and that from the g reconstruction 
decreases to the order of 10"^. We notice that with the noise 
from intrinsic ellipticities included, the correction due to the 
smoothing order problem gets larger. For our peak counts 
analyses here, we have the base shear and convergence maps. 
Therefore we can construct the direct noise-added shear and 
convergence fields and then calculate a at each grid point. 
With the a field, the correction Ax at each grid point can be 
obtained. The corrected results are shown in Figure 8. The 
blue histogram is for the corrected peak counts, which is very 
consistent with the results from either the direct noise-added 
convergence maps (upper left panel of Figure 2) or the shear 
reconstructed maps (left panel of Figure 5), and is in good 
agreement with the model prediction (red histogram). The 
green histogram is for the uncorrected peak counts. The right 
panel shows the parameter constraints from the corrected peak 
counts. The meaning of the symbols and lines are the same 
as that in the right panel of Figure 6. It is seen that the best 
fit values (red symbol) are very close to the best fit from the 
shear reconstruction (green symbol) and agree with the fidu- 
cial ones (blue symbol) very well. Thus for clarity, in the fol- 
lowing discussions of the mask effects, we use the corrected 
convergence maps from g reconstruction in both cases with 
and without masks. 

It should be pointed out that in real observational analyses, 
we cannot perform such corrections for the smoothing order 
problem straightforwardly because we do not know the true 
signals. We will explore ways to handle such effects in our 
future studies. 

4.2. Mask effects 

In this section, we discuss the mask effects on weak lens- 
ing peak analyses by comparing two sets of ^-reconstructed 
convergence maps corrected for the smoothing order prob- 
lem with and without masks, respectively. There are totally 
160 maps for each set. For each map in the case without 
masks, there is a corresponding map with masks that the 



source galaxies are exactly the same as the other one except 
that the galaxies within the masked regions are discarded. We 
then have 160 pairs of maps that allow us to do detailed com- 
parisons. The mask model is described in §3.3. Three cases 
with the number of masks N„,ask = 140, 280 and 420 for each 
9 deg^ (corresponding mask fraction of ~ 6%,~ 13% and 
~ 19%, respectively) are considered and Nmask = 280 is taken 
to be our fiducial case for most of the results presented in the 
following. 

To perform detailed studies for peaks in convergence maps 
with and without masks, we need to identify the peak corre- 
spondences between each pair of maps. This is done by peak 
matching. For each peak in a map from one set, we check 
for peaks within 3.5 arcmin in each dimension around it in its 
peer map from another set and define the nearest peak within 
this region as its partner peak. Only those pairs of peaks that 
are partners to each other are identified as peaks with corre- 
spondences. 

4.2. 1 . Effects on weak lensing peak analyses 

The existence of masks affects the convergence reconstruc- 
tion and consequently the peak properties both in their spatial 
location and peak height. 

In Figure 9 we present the mask effect on spatial positions 
of peaks for Nmask - 280. The left panel shows an exam- 
ple map of the spatial distribution of peaks with correspon- 
dences. The squared and plus symbols are respectively for 
peaks in the cases with and without masks. The red, black 
and yellow ones are for the pairs of peaks with their spatial 
offset larger than 0.5 arcmin, in the range of [0.2, 0.5] arcmin 
and less than 0.2 arcmin, respectively. It is seen clearly that 
strongly affected peaks are almost all closely associated with 
masks, especially large masks. In the right panel of Figure 
9, the statistical offset distribution averaged over 160 pairs of 
maps is shown. There are about 40% of peaks with offset 
larger than 0.1 arcmin. The fraction with offset larger than 
0.5 arcmin is ~ 11%. We also notice that lower peaks are 
more strongly affected by masks. This offset due to mask can 
have significant effects on weak lensing analyses for individ- 
ual clusters. For a typical weak lensing observation targeting 
at a particular cluster, the observed size is about 20 arcmin. If 
there happens to be a large mask close to the central region of 
the cluster, the weak lensing determined center for the clus- 
ter can be considerably offsetted from its true center, which in 
turn can lead to large errors in the weak lensing determination 
of the density profile for the cluster. 

We now discuss the mask effects on peak heights. Figure 10 
shows the results, where the left panels are the peak number 
distribution in 9 deg^ averaged over 160 maps in each case 
and the right panels are the peak number differences between 
the cases without and with masks. The upper and lower pan- 
els are for the peak height measured in K^ and in signal-to- 
noise ratio V - Kn/ctq with cro = 0.02, respectively. It is 
clearly seen that the number of peaks in high signal-to-noise 
bins is systematically higher in the case with masks, which 
can expectedly affect the cosmological parameter constraints 
with weak lensing peak counts significantly. We further exam 
the correlation between the positions of the strongly affected 
peaks and the locations of masks. We define two types of 
strongly affected peaks. Type I is for peaks with their peak 
height difference between the cases with and without masks 
higher than Icr. Type II is for peaks without correspondences 
between the two cases. Figure 11 presents a typical map with 
masks. The left panel shows the spatial distribution of peaks 
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with squares and pluses for peaks in the cases with and with- 
out masks, respectively. The red symbols are for Type 1 peaks, 
the green symbols are for Type 11 peaks with the purple ones 
for Type 11 peaks with v > 2, and the yellows ones are for 
the rest. The clustering of the strong affected peaks around 
large masks are apparent. The right panel shows particularly 
the Type 1 peaks with v < 3.25 in the case without masks but 
with the corresponding peak height shifting to v > 4.25 in the 
case with masks. It is found that they all trace large masks. 
It is these high Type 1 peaks that can affect profoundly the 
cosmological parameter constraints discussed here. 

Figure 12 shows the fitting results. The peak counts av- 
eraged over 160 maps in each case are shown in the upper 
left panel. The blue, green and red histograms are for peak 
counts in the case with masks, without masks and the theoret- 



ical prediction of Fan et al. (2010) with a uniform noise with 
ctq - 0.02. It is seen clearly that the peak counts considered 
here are systematically larger due to the presence of masks. 
The upper right panel shows the fitting results with the survey 
area of 9 deg^. The fitting is done with the 'observed data' 
being the peak counts for the masked case, and the model of 
Fan et al. (2010) with cr^) = 0.02. The meanings of the lines 
and symbols are similar to those of Figure 3. Clearly, the en- 
hanced peak counts due to the occurrence of masks lead to a 
large bias in cosmological parameter fitting. Even for a sur- 
vey of 9 deg^, the true cosmological parameter values (blue 
symbols) lie outside the 2cr contour around the best fit (red 
symbol). This demonstrates the significance of the mask ef- 
fects, which must be taken into account carefully in cosmo- 
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Figure 11. The illustration of the spatial distribution of the aifected peaks. In the left panel, the red, green and purple symbols are for Type I aifect peaks. Type 
II affected peaks and Type II with v > 2, respectively. The yellow symbols are for the rest of the peaks with correspondences between the case with and without 
masks. The right panel shows the Type I peaks with v < 3.25 in the case without masks but with the corresponding peak height shifting to v > 4.25 in the case 
with masks (red symbols). 



logical parameter constraints with weak lensing peak counts. 
The lower left and right panels show the results for the survey 
areas of 100 deg^ and 18000 deg^, respectively. The bias due 
to the mask effects is intolerable. 

The above results are shown for the average masked area 
fraction of ~ 13% with N„,ask - 280 in 9 deg^. We also ana- 
lyze how the effects depend on the masked fraction. We con- 
sider three cases with the number of masks Nmask - 140, 280 
and 420 in 9 deg^ and the corresponding masked fraction of 
~ 6%, ~ 13% and ~ 19%, respectively. The peak statistics 
are listed in Table 1 . The mask effects are clearly stronger for 
larger masked fraction. The fraction of peaks with v > 3 is 
about 9% in the case without masks. This fraction increases 
to ~ 11%, ~ 13% and ~ 15.7% for N,„a,k = 140,280 and 420, 
respectively. More than 90% and 70% of Type 1 and Type 
11 affected peaks, respectively, are within the regions around 
masks with a size of twice the mask radius. The results are 
further visually illustrated in Figure 13 with all the symbols 
the same as those shown in left panel of Figure 1 1 . The corre- 
sponding fitting results for the survey area of 9 deg^ are shown 
in Figure 14. We see that with the increase of the masked 
fraction, the effects become larger. For the masked fraction of 
~ 19%, the bias for (Q,„,cr8) is already larger than 3cr for a 
9 deg^ survey. 

4.2.2. Mask effects correction 

We have demonstrated in §4.2 that the mask effects on weak 
lensing peak counts are significant. The subsequent cosmo- 
logical parameter constraints are largely biased if they are not 
taken into account properly. We therefore need to explore 
ways for how to control the mask effects on cosmological ap- 
plications with weak lensing peak accounts. 

From Table 1 and Figure 11, we see that the strongly af- 
fected peaks are mostly clustered around masks. Thus the 
first method we use to suppress the mask effects is to exclude 
the severely affected regions around masks when preforming 
the peak counting. It is expected that the bias on cosmologi- 
cal parameters can be considerably removed but inevitably at 



the expense of loosing effective survey areas and therefore en- 
larging the statistical error contours. We name this method as 
the rejection method. We consider three cases with the rejec- 
tion regions of 1 , 1 .5 and 2 times of the mask size around each 
mask. We pay attention to the mask overlaps. Figure 15 shows 
the results. The upper left panel shows the peak number dis- 
tributions. The blue histogram is for the peak counts without 
applying rejection. The black one is for the peak counts after 
rejecting regions of 1.5 times of the mask size. The red one 
is the model prediction of Fan et al. (2010) with a uniform 
noise level of ctq = 0.02. All the peak counts are scaled to 
1 deg^. The upper right, lower left and right show the fitting 
results for the three considered rejections, from the smallest 
to the largest rejections, respectively. It is seen that rejecting 
only the masked areas is not enough (upper left panel) and 
the bias is still apparent. With the rejection of regions of 1.5 
times of the mask size around masks, the bias is suppressed 
to an insignificant level (lower left). To increase the rejection 
areas further does not improve the bias any more (lower right 
panel). It is also seen that with the increase of the rejection 
areas, the confidence contours become larger as expected. We 
thus conclude that to reject regions of 1.5 times of mask size 
around masks is an optimal choice that can remove the bias in 
cosmological parameter constraints without loosing statistics 
significantly. 

We now consider another approach in which we still keep 
the regions around masks in the peak counting. Here how- 
ever, the regions around masks and the regions away from 
masks must be theoretically considered differently. In other 
words, we need to revise our theoretical model. The major 
difference between the two types of regions is the effective 
number of galaxies that can be used in the convergence re- 
construction. In regions around masks, much fewer galaxies 
are usable. Therefore the noise level in those regions is sys- 
tematically higher. This higher noise affects the peak counts 
in regions around masks in two ways. One is that the system- 
atic peak hei ght shift for true p eaks is larger than regions ways 
from masks ( Fan et al. |2010[ l. The other is the enhancement 
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Figure 12. The bias in cosmological parameter constraints due to the mask effects. The upper left panel shows the peak number distribution where the blue, 
green and red histograms are for the cases with masks, without masks and the model prediction of Fan et al. (2010) with o"o = 0.02. The upper right panel 



shows the fitting results for the survey area of 9 
case with masks as the 'observed data' and ctq 
18000 deg , respectively. 



deg . The blue symbol is for the fiducial values and the red symbol is for the best fit with the peak counts in the 
= 0.02 in the model of Fan et al. (2010). The lower left and right panels are for the survey area of 100 deg and 



of the number of noise peaks given their peak heights mea- 
sured in signal-to-noise ratio with ctq still be 0.02, the value in 
regions away from the masks. The later can be understood by 
noting that for any two-dimensional Gaussian random fields, 
their peak number density distributions are the same if the 
peak height is measured in signal-to-noise rati o with the noise 
level (To different for different noise field (e.g., | van Waerbeke] 
2000 1. Here for regions around masks, their true noise level is 
high. When counting peaks, however, we still use ctq = 0.02 
(corresponding to rig - 30 arcmin"^) uniformly to define the 
signal-to-noise ratio of peaks. Therefore the number density 
of peaks with high v = Kn/ctq (cto = 0.02) is higher than 
that in regions away from masks. Based on these considera- 
tions, we develop a two-noise-level model that takes into ac- 
count the lower number density of usable galaxies in regions 
around masks. Specifically, for each grid point, we calculate 
the effective contribution of different galaxies to it by 



R' = 



YjI^1\ KMjt 



Ruk 






(27) 



where R - \ for galaxies outside masks and R - for galax- 
ies inside masks. The kernel ma is the Gaussian smoothing 
function with 6g = I arcmin, consistent with our analyses. 
We then obtain an effective number density of source galaxies 
at each grid point by 



nl = R'ng. 



(28) 



where n^ - 30 arcmin 



here. 



We divide the considered survey area into two parts, the 
regions around masks with 1.5 times of the mask size and 
the rest of the regions away from masks. They are named 
as, correspondingly, the mask affected regions and the nor- 
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Figure 13. The illustration of the dependence of the mask effects on the masked fraction. From left to right, the masked area fraction is ~ 6%, ~ 13%, and 
~ 19% (Wjiiflji = 140, 280, and 420 in 9 deg), respectively. The meanings of the symbols are the same as those in the left panel of Figure 11. 

Table 1 

Mask effects on peak statistics with different mask fractions. 



Jmasii ^masl Jnocon Joffsel /v>: Jm,v >1 JLPii ^ JNCii JLlinALL 



13% 
19% 



140 

280 
420 



6.07% 
10.76% 
15.02% 



24.56% 
40.07% 

52.25% 



8.83% 
).01% 
!.93% 



10.88% 
13.21% 
15.76% 



93.30% 
93.50% 
93.36% 



71.52% 41.02% 
77.39% 44.40% 
82.60% 47.90% 



"Masked area fraction 
''Number of masks in 9 deg 

"^Fraction of no-correspondence peaks among the total number of peaks 
Fraction of peaks with spatial offset larger than 0. 1 arcmin. 
■^Fraction of peaks with S/N > 3 in the mask-free case 
'Fraction of peaks with S/N > 3 in the case with masks 

^Fraction of Type I affected peaks within regions around masks with a size of twice the correspinding masks among the total number of Type I affected peaks 
''Fraction of no-correspondence peaks within regions of twice the size of masks among the total number of no-correspondence peaks 
'Fraction of (Type I+Type II) affected peaks within regions of twice the size of masks among the total number of peaks within the regions 



Mask fraction - 6% 



Mask fraction ~ 1 3% 



Maskfraction ~ 1 9% 
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Figure 14. The corresponding results of cosmological constraints for different masked fractions shown in Figure 13. 



mal regions. We calculate the effective number density of 
usable galaxies for the mask affected regions n^™"**^ by av- 
eraging n'L of all the grid points within these regions. We 
obtain n™""'' ^ 13.7 arcmin"^ for the case with the masked 
fraction of 13% considered here. For the normal regions, the 
effective number density n"'""" is calculated similarly by av- 
eraging all the grid points in the normal regions. We have 
n™""" * 29.3 arcmin"^, very close to «g = 30 arcmin"^ as 
expected. With these two effective number densities of galax- 
ies, we calculate the number density of peaks separately for 
mask affected regions and for normal regions with the model 
of Fan et al. (2010). We then rescale the height for peaks 
by a uniform noise level cro - 0.02, the value used in count- 
ing peaks from simulation maps. The theoretical predictions 
for the number of peaks in each S/N bin (with ctq = 0.02) 



in the mask affected regions and normal regions are calcu- 
lated by multiplying the corresponding areas of the two re- 
gions with their number densities of peaks. We then sum up 
the peak numbers in each bin in the two regions to obtain 
the total number of peaks in each bin. These theoretical pre- 
dictions are then compared with the 'observed' numbers of 
peaks directly from the convergence maps with masks with- 
out any area exclusions. The results are shown in Figure 16 
with the left panel for the peak counts and the right panel for 
the fitting results from our two-noise-level model. The blue, 
black and red histograms in the left panel correspond to the re- 
sults from convergence maps in the case with masks directly, 
the theoretical prediction from the two-noise-level model, and 
the theoretical prediction from the model of Fan et al. (2010) 
with a uniform noise level of ctq = 0.02. Comparing to the 
upper right panel of Figure 12, the fitting results improve dra- 
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Figure 15. Results after rejection of regions around masks. Tlie upper left panel shows the peak number distribution per deg with the blue, black and red 
histograms for peak counts without rejection, with rejection of regions of 1.5 times the mask size around masks and the model prediction of Fan et al. (2010). 
The upper right panel shows the fitting results with the rejected regions equal to the masked regions. The lower left and right panels are the results for rejections 
of regions with 1.5 and 2 times of mask size, respectively. Here the number of masks is Npg„;; = 280 in 9 deg and the corresponding masked fraction is 13%. 



matically with a much reduced bias considering the parameter 
degeneracy along the dotted line. The best fit;^'^ » 0.16 here 
while it \s x^ ~ 6.64 with the single noise model. The results 
show that our two-noise-model works very well and can po- 
tentially be applied in cosmological studies with weak lensing 
peak counts in the presence of masks. 

5. SUMMARY AND DISCUSSION 

In this paper, we analyze the mask effects on weak lensing 
convergence peak statistics and the consequent cosmological 
parameter constraints from weak lensing peak counts. Along 
the way, we also study different effects arising from conver- 
gence reconstructions from measured ellipticities of source 
gala xies. Our base conve rgence and shear maps are taken 
from [White & Vale | ( |2004j with the source redshift z., = 1. 
We generate source galaxy catalogs by randomly populating 
galaxies on the source plane with rig - 30 arcmin"^. A ran- 
dom intrinsic ellipticity based on the truncated Gaussian dis- 



tribution with the total cr^ - 0.4 is assigned to each galaxy. 
The mask catalogs are produced according to the mask size 
distribution of CFHTLS. To perform statistical analyses, for 
each of the 32 base weak lensing simulations, we generate 5 
noisy maps without/with masks with different assigned intrin- 
sic ellipticities and different masks in the masked case. We 
therefore have totally 160 convergence maps for each case. 
We use the theoretical model of Fan et al. (2010) for weak 
lensing peak counts to analyze the cosmological parameter 
constraints in different cases. The model takes into account 
fully the noise effects arising from intrinsic ellipticities so that 
we can use directly peaks from convergence maps in cosmo- 
logical studies without the need to find true peaks associated 
with clusters of galaxies. We test the model applicability with 
the 'observed data' from simulations without masks and find 
that the derived cosmological parameters agree with the fidu- 
cial input of the simulations very well without significant bias. 
Our main results are summarized as follows. 
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Figure 16. Left: Cosmological constraints for ~ 13% masked g reconstruction after equivalent correction. Same notation as upper-right panel of Fig 2. Right: 
Average peak height distribution for masked g reconstruction (blue), equivalent corrected model (black) and fiducial model prediction (red). 



(1) Because of the difference between the smoothed (g) and 
(t)/(1 - {k}), the smoothing procedure in g reconstruction 
leads to the reconstructed convergence systematically lower 
than the true smoothed convergence in dark matter halo re- 
gions. This can affect the weak lensing peak counts, which 
should be considered carefully, especially for future studies 
with large surveys. 

(2) The mass-sheet degeneracy problem in convergence re- 
constructions can affect the covariance matrix of peak counts, 
which in turn affects the confidence range of cosmological pa- 
rameter constraints. It is expected that with the increase of the 
area in reconstruction, this effect in covariance matrix should 
decrease. 

(3) The occurrence of masked regions reduces the number 
of usable source galaxies and therefore increases the noise 
in the regions around masks. This in turn leads to system- 
atic increases of the number of high peaks and consequently a 
significant bias in cosmological parameters constrained from 
weak lensing peak counts. The larger the masked area fraction 
is, the larger the effects are. 

(4) We find that the significantly affected region around a 
mask is about 1.5 times of the mask size. Excluding such re- 
gions in peak counting can eliminate largely the mask effects 
and reduce the bias in cosmological parameters derived from 
weak lensing peak abundance significantly. 

(5) We develop a two-noise-level model that treats the mask 
affected regions separately. This model can account for the 
mask effects on weak lensing peak counts very well. The de- 
rived cosmological parameters based on this model improve 
dramatically comparing to the large bias from the model with 
a uniform noise. 

In our analyses, we apply the Kaiser-Squires method with 
the nonlinear reduced shear in the convergence reconstruc- 
tion. For the maximum-likelihood reconstruction method 



( Bartelmann et al. ||1996| , we expect that the mask ef- 
fects on the reconstructed convergence field and the peak 
counts are qualitatively similar to the results shown in this 
paper although quantitative studies are still needed. For other 
methods, such as the multi-scale entr opy restoration filtering, 
namely MRLens ( Starck et al. 2006|l, the mask effects can be 



different and detailed analyses should be carried out when a 
specific reconstruction method is used. 

We note that the two-noise-level model works well to take 
into account the mask effects on weak lensing peak counts. 
Detailed inspection finds, however, that the model predic- 
tion for the dependence of peak counts on the peak height 
is slightly steeper than that from simulations. To improve the 
model, we may need to treat the mask affected regions more 
accurately by further dividing them into different sub-regions 
with different noise levels. We will investigate along this line 
of approach thoroughly in the future. 
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